Strain Gradients in Epitaxial Ferroelectrics 
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X-ray analysis of ferroelectric thin layers of Bai/2Sri/2Ti03 with different thickness reveals the 
presence of internal strain gradients across the film thickness and allows us to propose a functional 
form for the internal strain profile. We use this to calculate the direct influence of strain gradient, 
through flexoelectric coupling, on the degradation of the ferroelectric properties of thin films with 
decreasing thickness, in excellent agreement with the observed behaviour. This work highlights the 
link between strain relaxation and strain gradients in epitaxial films, and shows the pressing need 
to avoid strain gradients in order to obtain thin ferroelectrics with bulk-like properties. 

PACS numbers: 77.55.-(-f, 77.80.-e, 68.55.-a, 61.10.-i 



The interest in ferroelectric thin films is rapidly ex- 
panding due to the recent development of both experi- 
mental techniques and calculation tools that allow to ex- 
plore the ferroelectric phenomena at atomic level 0, . 
The incorporation of realistic mechanical and electrical 
boundary conditions in the first-principles formulations 
is generating new insight on the mechanisms limiting the 
ferroelectric response in thin ferroelectric layers H, ^ |^ ■ 
But while the new evidence suggests that ferroelectricity 
may indeed be stable in thin films only a few monolay- 
ers thick P, |3j , the sharp peak in dielectric constant 
usually associated with the ferroelectric transition is sys- 
tematically depressed in thin films. This obviously limits 
the technological impact that would arise from the ability 
to maintain ferroelectricity and large dielectric constants 
down to the nanoscale in real devices. 

Strain caused by lattice mismatch with the substrate is 
an important factor affecting the properties of thin films. 
Strain can modify the phase diagram of epitaxial ferro- 
electrics 0, 01, change the order of the transition^ 
and shift transition temperatures 0)0- However, strain 
alone does not generally account for the observed smear- 
ing of the dielectric peak, as a sharp anomaly is still 
expected at the strain-modified transition temperature. 
Gradient terms (of strain, composition, defects, etc) have 
recently been propos ed to account for the reduced dielec- 
tric constant However, no experimental studies 
have provided quantitative insight in the gradient terms. 
The aim of the present work is to detect and measure 
strain gradients in a set of lattice-mismatched epitax- 
ial thin films, and to correlate the measured gradients 
with the measured dielectric properties. The tools used 
in this work can be applied to any material where gradi- 
ents, not just of strain but also of impurities or vacancies, 
are expected play a role. Showing the link between strain 
relaxation and strain gradients has therefore wider impli- 
cations beyond ferroelectricity and is an important result 
for general thin film epitaxy. 




FIG. 1: (Color) Williamson-Hall plots as a function of film- 
thickness. Inset: Diffraction peaks show that broadening is 
maximum at intermediate thickness. The peaks at 46.1 deg. 
correspond to the electrode (SRO) 



The films studied in this work are Bao.sSro.sTiOa 
(BST) dielectric layers with SrRuOs (SRO) bottom elec- 
trodes, grown by pulsed laser deposition. Epitaxy has 
been verified by cross-sectional high-resolution transmis- 
sion electron microscopy (TEM). Details of the growth 
and TEM characterisation are published elsewhere |3| . In 
the present work the crystallographic analysis has been 
performed using a Philips X'pert MRD diffractometer 
with CuKairadiation (A=1.540 A). 

The lattice parameters are extracted from the posi- 
tion of the perovskite pseudocubic (002) diffraction peak 
(see inset of Figure 1). This allows the calculation of 
the average out-of-plane strain in each film, given by 



_ c{t)-co 



, where c is the average out-of-plane lattice 
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parameter, Cq is the reference value, i.e., the bulk lattice 
parameter {cbst= 3.95 A), and t is the film thickness. 
The average strain for each film is shown in the inset 
of Figure 2. There is an out-of-plane expansion for the 
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thinnest films that decreases with thickness, as observed 
before |3, ■ This is consistent with relaxation of the 
in-plane compression induced by the smaller lattice pa- 
rameter of the bottom electrode {csro =3.93A). 

Similar to what is known for semiconductor and metal- 
lic epitaxial layers, perovskite oxides are known to relieve 
strain as film thickness is increased. The strain-relieving 
mechanism is thought to be mainly the formation of mis- 
fit dislocations. As thickness increases, the accumulation 
of elastic strain energy overcomes the barrier for the for- 
mation of misfit dislocations, which ease the strain. It is 
generally implied that this strain relaxation is homoge- 
neous across the film, and thus the strain state is only a 
function of the film's total thickness: e=e(t). 

However, around a dislocation the lattice is locally dis- 
torted. The accumulation of misfit dislocations at the 
film-substrate interface means that the film's lattice pa- 
rameters are more distorted near the clamped surface 
than at the free one, leading to stress/strain distributions 
0,0,0,0. Furthermore, strain may not be relaxed 
solely by dislocations, as other inhomogeneous mech- 
anisms (such as vertical segregation of different-sized 
cations) have been observed Thus, rather than a 

quantity dependent only on the thickness i, strain should 
be described as an internal profile dependent also on the 
distance to the film-substrate interface, z: e=e{z,t). 

The homogeneous vs inhomogeneous scenarios of strain 
relaxation are not only different from a structural point 
of view, but have consequences for the functional prop- 
erties. Inhomogeneous strain fields around dislocations 
[TsL 01 and impurities 0| affect the polarisation and 
critical temperatures of ferroelectric thin films. Crucially, 
also, inhomogeneous strain is necessarily associated with 
local strain gradients, which couple to the polarisation 
via the flexoelectric effect 0,0,120^. Measuring the ver- 
tical strain gradient is therefore essential to correctly de- 
scribe the functional properties of ferroelectric thin films. 

In order to calculate the strain gradients, x-ray diffrac- 
tion peak broadening has been analysed as a function 
of film thickness. There are at least two contributions 
to peak broadening: one due to the finite thickness 
of the sample, and another due to the inhomogeneous 
strain. The two have different angular dependence, and 
can therefore be separated by looking at peak width 
for different reflections and fitting the results using the 
Williamson-Hall relation l21l: 
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(3cose^ K^+Ae, sin(6'). 



(1) 



where D is the coherent length perpendicular to the film's 
surface (roughly proportional to the film's thickness), A 
is the X-ray wavelength (A=1.54A in our case), 6 is the 
diffraction angle, (3 is the peak integral breadth (close 
to the full width half maximum) minus the instrumental 
broadening, and K is an empirical constant close to 1. 

Linear fits of (3cos9 vs sin6' yield the coherent length D 
and inhomogeneous strain for each film. We have per- 
formed such fits for the (OOh) (h=l:4) reflections, find- 
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FIG. 2: Inhomogeneous strain as a function of film thickness. 
The solid line is a visual guide. Inset: average out-of-plane 
strain. Dotted lines are least-square fits using the averages of 
Eqs. Hand El 



ing the linear dependence excellent for all our samples 
(r^>0.9)(Figure 1). It is nevertheless worth mentioning 
that the linear W-H plots are only one of existing strate- 
gies to separate size and strain broadening. Quantitative 
results for the inhomogeneous strain may therefore vary 
depending on the approach used 22]. 

Figure 2 shows the inhomogeneous strain as a func- 
tion of film thickness. The existence of a maximum in 
ei confirms the presence of a monotonically decreasing 
internal strain profile as a function of z in the films. For 
very thin films the reduced thickness means a small dis- 
persion in lattice parameters; conversely, for very thick 
films there may be a large difference between the lat- 
tice parameters at the clamped and free interfaces, but 
the volume fraction of totally relaxed material is large 
and dominates the diffraction, so that again the variance 
is small. In between, there is an intermediate thickness 
where there is a large dispersion in lattice parameters, 
with their thickness fractions being similar. At this point 
the inhomogeneous strain is maximum. 

Establishing the actual mechanism of strain relaxation 
is beyond the scope of this paper, but it is worth men- 
tioning that the observed behaviour is consistent with 
predictions for dislocation-based relaxation: the diffrac- 
tion peak width in this case is expected to be propor- 
tional to y/p/t, where p is the linear dislocation density, 
which grows rapidly around a critical thickness and then 
saturates psj . 

Extracting quantitative values for the internal strain 
profile from this analysis requires solving the integral 
equations for the average (e) and inhomogeneous (e^) 
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strain: 
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where e, are extracted from peak position and width, 
respectively, and e(z, t) = <:{^^*-)~'^o ^-^^ internal strain 
profile. 

The easiest way to resolve these equations is the inverse 
approach: assume a certain shape for the internal strain, 
solve the integrals JSJ, ®, and modify the functional pa- 
rameters to achieve a good match with the experimental 
results. This method relies heavily on the correct choice 
of functional dependence for e(z, t). As such, the re- 
sults of the quantitative analysis should be treated only 
as approximations. 

A general model for the strain profile, independent of 
the actual relaxation mechanism, reflects that strain re- 
laxation should be proportional to the strain itself, which 
yields an exponential dependence on z |24| : 



9e 
dz 



e(z) = eoe 



(4) 



where eg is the strain at the film-substrate interface, and 
(5 is a measure of the penetration depth of the strain. If 
dislocations are considered as the main relaxation mech- 
anism, a recent strain-gradient theory predicts the 
vertical profile in the layers to be given by 



e(z, t) = eo 
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It is worth noticing that Eq. Q) is a limiting case 
of lO when the film thickness is larger than the strain 
penetration depth (t:^ 6). 



Either of these expressions can be integrated to yield 
e{t). The least-squares fits to the experimental results 
using both are shown as dashed lines in the inset of figure 
1. The value of the fitting parameters is eo = 0.013 ± 
0.001 and (5 = 60 ± 12 nm, and eo=0.010±0.001 and 
(5=85±13 nm for the fits with the averages of Q and © 
respectively. 

The calculated curves for e^, using the above parame- 
ters also reproduce the experimental results for the max- 
imum value and associated thickness. However, beyond 
the maximum, the predicted relaxation of Ci using these 
equations is slower than experimentally measured. This 
discrepancy may be explained by the presence of more 
than one relaxation mechanism, each with different pen- 
etration length S. Furthermore, eg is a function of t, 
since the increase in dislocation density for thicker films 
affects the strain at the film-substrate interface. Thus, 
both parameters should, in principle, be considered as 
thickness-dependent: eo{t), S{t). In order to calculate 
the thickness dependence of eo(^) and S{t) we note that 
there are two parameters and two equations to describe 
e and e^, so it is possible to calculate eo and S for each 
film separately. We have done this for the exponential 
strain profile Combining the Eqs. (0) and (PJ we 
can eliminate eo: 
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This is solved for each film in order to find 5{t), which 
is then used to calculate eQ(t). 

Once the internal strain profile e(z,t) is known, the 
strain gradient contribution to the functional properties 
can be calculated using an elastodielectric free energy 
expansion incorporating the fiexoelectric contribution: 
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where P is the out-of-plane polarisation; Sy the elastic 
compliances; a the in-plane stress (related to the mea- 
sured out-of-plane strain by the Young's modulus and 
Poisson's ratio: a = eY{-2v))] Qis is the transverse elec- 
trostrictive coefficient, C and D are the constants related 
to the energy contributions from polarisation and stress 
gradient, and 7 and u are, respectively, the direct and 
converse fiexoelectric coefficients. P is calculated by min- 
imising the thermodynamic potential, while the second 
derivative of G with respect to P yields the inverse per- 
mittivity. This is averaged over the thickness of the film 



to yield the effective value. The values of the coefficients 
used in this expansion are the same as in 10] . 

The relative dielectric constants calculated using the 
strain gradient extracted from our crystallographic anal- 
ysis are shown in Fig. (2), along with experimentally mea- 
sured for the same set of films. The predicted and mea- 
sured temperatures of maximum permittivity (Tm) are 
shown in the inset. 

Clearly, the decrease in dielectric constant and up- 
ward shift of T„j are well reproduced. Quantitatively, 
the prediction for T^a as a function of thickness is re- 
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FIG. 3: Calculated and measured relative dielectric constant 
as a function of temperature for films of thickness 660 nm, 
340nm, 280nm, 220nm and 145nm. Inset: Temperature of 
maximum permittivity, experimental (dots) and calculated 
(solid line) 

markably good, while the calculated dielectric constant 
is larger than experimentally measured. This was ex- 
pected, as our model does not take into consideration any 
other permittivity-depressing factors. The results show 



that the contribution to the depression in permittivity 
with decreasing thickness (size effect) from flexoelectric- 
ity alone is enormous. This is particularly valid when 
compared with the huge permittivities recently measured 
in gradient- free ferroelectric films 25] . 

It is worth noting that the dielectric constant is low- 
est for the thinnest films in spite of the relatively small 
value of ti . This is a natural consequence of the fact that 
the size effect is not caused by the inhomogeneous strain 
itself, but by the strain gradient, which is largest for the 
thinnest films. We emphasise also that while compressive 
in-plane strain can indeed be used to stabilise the ferro- 
electric state, this may come at the expense of reducing 
the permittivity if strain gradients are not avoided. Also, 
the procedure outlined here could be used to estimate 
vertical oxygen vacancy distributions, or gradients due 
to impurity concentration. Finally, these methods open 
the scope for studying the effect of strain gradients on 
other functional materials. 

In summary, X-ray analysis of peak broadening as a 
function of thickness shows that relaxation of strain in 
epitaxial films is associated with the appearance of inter- 
nal strain gradients. The dielectric constants calculated 
using these strain gradients are close to experimentally 
measured, clearly showing the fundamental role plaid 
by flexoelectric coupling in decreasing the dielectric con- 
stant. This work shows the urgent need to avoid strain 
gradients in order to prevent degradation of the ferro- 
electric response in thin films. 

Useful discussions with T. Hibma, E. Van der Giessen, 
T. Palstra and D. Boer are gratefully acknowledged. 



R. Ramesh and D.G. Schlom, Science 296, 1975 (2002). [14 
C.H. Ahn, K.M. Rabe and J-M Triscone Science 303, 
489 (2004) and refs therein. [15 
J. Junquera, R Ghosez, Nature 422, 506 (2003). 
I. Kornev, H. Fu, L. Bellaiche, Rhys. Rev. Lett. 93, [16 
196104 (2004). \\1 
C. Bungaro, K.M. Rabe, arXiv:cond-mat/0410375j 
O. Dieguez et al. Rhys. Rev. B 69, 212101 (2004). [IS; 
N. A. Pertsev, A. G. Zembilgotov, A. K. Tagantsev, Rhys. [19 
Rev. Lett. 80, 1988 (1998) 
C. Basceri, S. K. Streiffer, A. I. Kingon, R. Wasser, J. [20 
Appl. Rhys. 82, 2497 (1997) [21 
L. Sinnamon, R.M. Bowman, J.M. Gregg Appl. Rhys. [22 
Lett. 81, 889 (2002). 

G. Catalan, L. Sinamon, M. Gregg, J. Rhys. Cond. Mat. [23 
16, 2253 (2004). 

A. M. Bratkovsky, A. R. Levanyuk, arXiv: [24 
cond- mat/0402100 (2004). 

C.L.Canedy et al. Appl. Rhys. Lett. 77, 1695 (2000). [25 
M.-W. Chu et al. Nature Materials 3, 87 (2004). 



S. R. Alpay, I. B. Misirlioglu, V. Nagarajan, R. Ramesh, 
Appl. Rhys. Lett. 85, 2044 (2004). 

L. Nicola, E. Van der Giessen, M. E. Gurtin, submitted 
to J. Mech. Rhys. Solids. 

J.-L. Maurice et al, Rhil. Mag. 83, 3201 (2003). 

D. Balzar, R. A. Ramakrishnan, A. M. Hermann, Rhys. 
Rev. B 70, 92103 (2004). 

Sh. M. Kogan, Sov. Rhys. Solid State, 5 2069 (1964) 
W. Ma, L.E. Cross Appl. Rhys. Lett. 79; 4420 (2001); 
ihid 81, 3440 (2002) 

A. K. Tagantsev Rhys. Rev. B 34, 5883 (1986). 

G. K. Wilhamson, W. H. Hall, Acta Metall. 1, 22 (1953). 
J. G. M. Van Berkum, R. Delhez, Th. H. De Keijser and 

E. J. Mittemeijer, Acta Cryst. A 52, 730 (1996). 

V. M. Kaganer, R. Kohler, M. Schmidbauer, R. Opitz, 

B. Jenichen, Rhys. Rev. B 55, 1793 (1997). 

H. Joon Kim, S. Hoon Oh, Hyun M. Hang, Appl. Rhys. 
Lett. 75, 3195 (1999). 

M. M. Saad et al, J. Rhys.; Cond. Mat. 16, 451 (2004). 



